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^^ \ We present a combined analytical and numerical approach based on the 

[^^ ■ Mori projection operator formalism and Monte Carlo simulations to study 

0^ . surface diffusion within the lattice-gas model. In the present theory, the 

average jump rate and the susceptibility factor appearing are evaluated 
through Monte Carlo simulations, while the memory functions are ap- 
proximated by the known results for a Langmuir gas model. This leads 
^2 \ to a dynamical mean field theory (DMF) for collective diffusion, while 

Q I approximate correlation effects beyond DMF are included for tracer dif- 

O ■ fusion. We apply our formalism to three very different strongly interacting 

systems and compare the results of the new approach with those of usual 
Monte Carlo simulations. We find that the combined approach works very 



^ ' well for collective diffusion, whereas for tracer diffusion the influence of 



interactions on the memory effects is more prominent. 
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1 Introduction 



Most of the existing tlieoretical investigations of surface diffusion ^^^ are 
for cases where the interparticle interactions do not play an important role. 
Approximations for strongly interacting systems have only been developed 
for some special cases PHTHj , but a general understanding of the interaction 



effects is still lacking. At finite coverages, there are two different diffusion con- 
stants. The tracer diffusion coefficient Dt is directly related to the motion 
of a tagged particle as observed in STM and field ion measurements, while 
the collective diffusion coefficient Dc describes the macroscopic density fluc- 



tuations as measured in field emission and optical grating experiments |[T9|| . 



Already in 1981, Reed and Ehrlich [10| tried to relate Dc also to the Individ 



ual adparticle motion. They proposed that Dc can be expressed as a product 
of a thermodynamical factor and an effective jump rate of adparticles. How- 
ever, there was no rigorous theoretical basis for this decomposition, in which 
the effect of dynamical correlations (memory effects) is neglected PDI . These 



memory effects originate from the fact that a backward jump of a diffusing 
particle is more likely than jumps to other sites. 

The purpose of this Letter is to present a new method for studying surface 
diffusion within the lattice gas model |2I]. It is based on a combination of an 



analytic approach and Monte Carlo (MC) simulations. The only approxima- 
tion involved in the new approach is that the memory effects are approximated 



by the known results for the Langmuir gas model |2^. In the Langmuir gas, 
there are no direct interactions, but the double occupancy of lattice sites is 
excluded. It turns out that our result for Dc is exactly of the form proposed in 
Ref. [|T^. We apply the theory to three different physical systems: O/W(110) 
], adlayer on a stepped substrate P^, and flexible, chainlike molecules on 



a flat substrate [|18[. We flnd that the predictions of the theory for collec- 



tive diffusion compare very well with full MC simulations, thus justifying the 



validity of the Reed-Ehrlich description |]T0[. Moreover, the new method is 
computationally very efficient. For tracer diffusion, however, we flnd the effect 
of interactions on the memory effects to be more pronounced. 



2 Theoretical description of surface diffusion 

As a starting point for studying collective diffusion, we focus on the density- 
density autocorrelation function S'(|r—r |,t) = {5n{r,t)5n{Y , 0)), where 5n(r,t) 
n(r,t) — {n{r,t)) with an occupation variable n{r,t) = 0, 1 at a lattice site r 
at time t. The corresponding Laplace- Fourier transform is denoted by ^(q, z). 
Using the Mori projection operator formalism [^, it can be shown that [jl5 



S(q, z) 



^-Kq)x(q)~' + M(q,z) 



where 6(q) contains microscopic jump rate information, x(q) is the thermo- 
dynamical susceptibility, and the memory function M(q, z) contains the dy- 



nainical correlations. The collective diffusion constant can be obtained from 
the correlation function S'(q, z) by examining its pole in the limits q — *> and 
z ^ [ p!^ , p!5| ] . In the limit q —^ 0, 6(q) ~ NJ2a ^aQa^ay where N is the num- 



ber of adparticles in a system, F^ is the average jump rate and ia is the jump 
length along the direction a = x,y. Also Xo = liKig->oX(q)/^ = (('^^)^)/(^) 
is just the compressibility of the adsorbate overlayer. Note that when the 
memory function M is left out, this constitutes a dynamical mean field theory 
(DMF) for the collective diffusion constant. 

Our novel combined approach consists of evaluating the average jump rate F^ 



and the compressibility Xo ^ot through further analytic approximations |15], 
but rather by direct MC simulations. Fq, is directly obtained from the suc- 
cess ratio of individual particle jumps in the canonical ensemble, while xo is 
more conveniently evaluated within the grand canonical ensemble. This pro- 
cedure is easily implemented for arbitrary interaction strengths and transition 



algorithms |21| within the lattice gas model. For the memory function M, we 



follow Ferrando et al. [p!3| , pr5| and approximate it by the known expression for 
the Langmuir gas model. Thus the effect of direct interparticle interactions on 
M is not taken into account. In the Langmuir gas model the memory function 
for collective diffusion is exactly zero p2| , p6| , and our approach corresponds to 
the DMF approximation for Dc- Eq. (Q) then leads to an expression for Dc 
exactly of the form first proposed by Reed and Ehrlich: 

f F 
4 Xo 



For tracer diffusion, a similar analysis can be performed for the self-correlation 
function S's(|r— r \,t) = {6ns{r, t)6ns{r , 0)) with 6ns{r,t) = ns{r,t) — {ns{r,t)), 
where the tracer occupation number at r at time t, ns{r, t) = 6{r — R(t)), now 
refers to a tagged particle at R(t). In the expression analogous to Eq. (|I]), the 
jump factor 6(q) remains the same, x{^) is replaced by unity and the memory 
function Ms(q, z) differs from the corresponding one for collective diffusion. 
Again, we approximate the correlation effects resulting from Ms{q^,z) by the 
known expression for the Langmuir model. The resulting expression for Dt 
becomes [O 



DtT^ = ffiO)T^ , (3) 



where f{6) is a known correlation factor ^-^ that depends only on the cov- 
erage 6 and the geometry of the lattice. To summarize, our approach yields a 
DMF result for the collective diffusion constant Dc- For the tracer diffusion 
constant Dx, approximate correlation effects beyond the DMF are included. 



3 Results for model systems 

We now apply our formalism to study diffusion for three different systems 
based on lattice-gas models. The first system is O/W(110) [^,^. We study 
this system at a coverage of ^ = 0.45 over a wide range of temperatures. 
At low temperatures, the system is in the ordered p(2 x 1) phase, while for 



T > Tc it is disordered. Details of the model and MC simulations of the 
diffusion coefficients are in Ref. |2^. In Fig. |l](a), we show the results for 
collective diffusion in an Arrhenius plot. The agreement between the present 
DMF result and the direct MC data is remarkably good. Furthermore, DMF 
works rather well in the strongly interacting region, i.e. in the ordered phase 
and even close to Tc, the difference being always less than 15%. For tracer 
diffusion shown in Fig. 0(b), on the other hand, the discrepancy between the 
MC data and the DMF corrected by the approximate correlation factor is 
more significant, and becomes most prominent in the ordered phase where the 
interaction effects are important. 

The second system we consider is an adsorbate layer on a substrate with 
equally spaced straight steps. In addition to repulsive nearest-neighbor inter- 
actions between the adsorbates, this model includes an extra binding energy 
at step edges, an extra barrier for climbing over step edges, and enhanced dif- 
fusion along step edges. A detailed description of the model and the formalism 



p8| is given in Ref. [0]. In Fig. ^ we compare the theoretical results for col- 
lective and tracer diffusion with the direct MC simulations at a temperature 
comparable to the adsorbate interactions and the various substrate-induced 
activation barriers. For collective diffusion, the DMF gives again a very good 
description. For tracer diffusion, with the chosen parametrization |23], the 



component along the direction perpendicular to the steps is still described 
rather well by the approximate theory. For the component along the direction 
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parallel to step edges, the approximate theory fails aheady at relatively low 
coverages. This can be explained by the increased concentration at step edges 
due to the extra binding energy there, which leads to nontrivial memory effects 
when tracer particles have to make detours via terraces to pass each other. 

The last system we consider, polymer chain molecules adsorbed on a flat 
substrate, is the most complex one. The chains are modelled by the two- 



dimensional fiuctuating-bond model ||2^, in which the polymer segments oc- 



cupy single sites on a lattice and occupation of nearest and next nearest neigh- 
bor sites is excluded. There are no direct interactions between the polymers in 
the model. However, an effective entropic repulsive interaction exists between 



the chains. A more detailed description of the model can be found in Ref. |18|. 
In Fig. |^(a), we show a comparison for the diffusion constants calculated from 
both the present theory and the full MC simulation as a function of concen- 



tration 6 [^ . Even for this complex system, there is good agreement between 
the two approaches for the collective diffusion. Again, for the tracer diffusion 
whose results are presented in Fig. |^(b), the discrepancy is significant, and 
the results from the two methods differ even qualitatively. 



4 Discussion and conclusions 

The three model examples above clearly illustrate that in the case of collective 
diffusion, dynamical correlations between successive jumps of individual par- 
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tides cancel out to a large degree even for a strongly interacting system, and 
the dynamical mean field theory (DMF) works well. For tracer diffusion, the 
memory effects depend more strongly on the interaction between the adparti- 
cles and the approximate correlation effects from the Langmuir gas model do 
not provide an adequate description. 

The success of the DMF for Dc has provided us with a better understand- 
ing of collective diffusion as well as a practical tool for its evaluation. From 
the theoretical viewpoint, Dc in the DMF description is related precisely to 
an average jump rate of the individual adparticles, thus justifying the phe- 
nomenological assumption of Reed and Ehrlich. It provides a conceptual link 
between the individual particle jump rates, such as those measured in STM 
and field ion microscope experiments, and macroscopic density fluctuation 
measurements like field emission and optical grating studies. On the practical 
side, the calculation of the DMF description, including the average jump rate 
and the compressibility of the adparticles, is much easier to evaluate numer- 
ically compared with a direct simulation of Dc- Our preliminary results for 
the three model cases indicate that this speedup is typically of the order of 
25 — 100, the method being most efficient for cases where large system sizes 
or unit cells are needed. 

To conclude, our combined method provides a powerful tool for studying re- 
alistic models of complex systems. At present, we are employing the DMF 
formalism to study polymer chain molecules with direct interactions between 



them and surfaces with high concentration of imperfections such as steps, 
kinks and traps. 
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Fig. 1. Results for (a) Dq and (b) Dx as a typical Arrhenius plot in the O/W(110) 
system at ^ = 0.45. The results of conventional MC simulations along the two 
principal directions x and y are shown by open symbols, while the corresponding 
results of the approximate theory are presented by lines. The critical temperature 
of the order-disorder phase transition is denoted by Tc- The error bars of the MC 
results are roughly of the size of the symbols for Dc and very small for Dt- 



Fig. 2. Diffusion constants (a) Dc and (b) Dt as a function of concentration 
for diffusion in direction perpendicular to the steps (x direction) and parallel to 
the steps {y direction) in a model for a submonolayer of interacting adsorbates 
on a stepped substrate. The MC results are denoted by the open symbols, the 
approximate theory is shown by lines, and all data for Dyy has been scaled by a 
factor of 0.02. The error bars of the MC results are roughly of the size of the symbols 
for Dc and very small for Dt- 
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Fig. 3. Results for (a) Dq and (b) Dt as a function of the concentration 6 in the 
model polymer system studied. The MC results for diffusion coefficients are given 
by squares, while the corresponding approximate descriptions are shown by lines. 
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